REMARKS ON THE ILL-POSEDNESS 
OF THE PRANDTL EQUATION 



DAVID GERARD-VARET, TOAN NGUYEN 

Abstract. In the lines of the recent paper [4], we establish various ill-posedness results 
for the Prandtl equation. By considering perturbations of stationary shear flows, we show 
that for some linearizations of the Prandtl equation and some C°° initial data, local in 
time C°° solutions do not exist. At the nonlinear level, we prove that if a flow exists in 
the Sobolev setting, it cannot be Lipschitz continuous. Besides ill-posedness in time, we 
also establish some ill-posedness in space, that casts some light on the results obtained by 
Oleinik for monotonic data. 



1. Introduction 
Our concern in this paper is the famous Prandtl equation: 

dtu + ud x u + vdyU — dyU + d x P = f, y > 0, 



(1.1) 



d x u + d y v = 0, y > 0, 
u = v = 0, y = 0, 
lim^+oo u = U(t,x), 

that was proposed by Ludwig Prandtl [12] in 1904 as a model for fluids with low viscosity 
near a solid boundary. This model is obtained formally as a singular limit of the Navier- 
Stokes equations, in the limit of vanishing viscosity. In this asymptotics, y = is the 
boundary, x is a curvilinear coordinate along the boundary, whereas u = u(t, x, y) and 
v = v(t, x, y) are the tangential and normal components of the velocity in the so-called 
boundary layer. The pressure P = P(t, x) and tangential velocity U = U(t, x) are given: 
they describe the flow just outside the boundary layer, and satisfy the Bernoulli equation 

d t U + Ud x U + d x P = 0. 

Finally, the source term / = f(t,x,y) accounts for possible additional forcings. We refer 
to [7] for the formal asymptotic derivation and all necessary physical background. We shall 
restrict here to two settings: 

• the initial value problem (IV P): 

(1.2) (t, x, y) in [0, T) x T x E+, n| t=0 = u (x, y). 

• the boundary value problem (BVP): 

(1.3) (i, x, y) in T x [0, X) x R+, u\ x=0 = m(t, y). 
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Although (jl.ip is the cornerstone of the boundary layer theory, the well-posedness of 
the equation and its rigorous derivation from Navier-Stokes are far from being established. 
The reason is that the boundary layer undergoes many instabilities, the impact of which 
on the relevance of the Prandtl model is not clear. One popular instability is the so-called 
boundary layer separation, which is created by an adverse pressure gradient (d x P > 0) and 
a loss of monotonicity in y of the tangential velocity: see [Tj- 

Roughly, only two frameworks have led to positive mathematical results: 

• the analytic framework: under analyticity of the initial data and the Euler flow, 
Sammartino and Caflish |13| showed the well-posedness of the initial value problem, 
and successfully justified the asymptotics locally in time. See also [9]. Up to our 
knowledge it is the only setting in which the Prandtl model is fully justified. 

• the monotonic framework: under a main assumption of monotonicity in y of the 
initial data, Oleinik and Samokhin [IT] proved the local well-posedness of both the 
boundary value and initial value problems. The latter result was extended to be 
global in time by a work of Xin and Zhang [14j when / = and d x P < 0. 

We refer to the review article [3] for precise statements and ideas of proofs. We stress 
that in all the aforementioned works, either analyticity or monotonicity of the initial data 
is assumed. When such assumptions are no longer satisfied, instabilities develop, and the 
Prandtl model is unlikely to be valid, at least globally in space time. For instance, it 
was shown by E and Engquist [2] that C°° solutions of (jl.ip do not always exist globally 
in time. As regards the asymptotic derivation of Prandtl, some counterexamples due to 
Grenier [5] have shown that the asymptotics does not hold in the Sobolev space W 1,00 . 
Finally, negative results have culminated in the recent paper (5| by the first author and 
Dormy, that establishes some linear ill-posedness for the initial value problem in a Sobolev 
setting. We shall come back to this article in due course. Broadly, the authors consider 
the linearized Prandtl equation around a non-monotonic shear flow, and construct 0{k~°°) 
approximate solutions that grow like e^ 1 for high frequencies k in x. Our aim in this note 
is to ponder on this construction to establish further ill-posedness results for the Prandtl 
equation. 

Let us now present our results. We only treat the case of constant U > 0, and shall 
consider perturbations of some steady shear flow solutions: 



Note that only non-trivial source terms / yield non-trivial solutions of this form. But 
up to minor changes, our results adapt to some unsteady shear flows (u(t,y),0) satisfying 
dtu — d^u = 0. Thus, the important special case / = can be treated as well. 
The system satisfied by the perturbation (u,v) of (u s ,0) reads: 



(u,v) = (u s (y),0) 



u s (0) = 0, 



lim u s = U. 



dtu + u s d x u + u' v + ud x u + vd y u — d% 




y >o 
y >o 
y = o 



(1.4) 



plus the condition lim y _ 5 . +00 u = 0, that will be encoded in the functional spaces. 
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Our first result is related to the linearized version of this system, that is: 

dtu + u s d x u + u' s v — dyU = 0, y > 0, 

(1.5) < d x u + d y v = 0, y > 0, 

u = v = 0, y = 0. 

We state a strong ill-posedness result for the initial value problem, namely: 

Theorem 1.1 (Non-existence of solutions for linearized Prandtl). There exists a shear flow 
u s with u s — U£ C^°(M_|_) such that : for allT > 0, there exists an initial data uo satisfying 

i) e y u G H°°(T x M+). 

ii) The I VP (jl.5[) - f)1.2|) has no distributional solution u with 

u G L°°(]0,r[;L 2 (T x M + )), d y u G L 2 (]0,T[xT x R+). 
We quote that a solution u of ()1.5p with the above regularity satisfies 

a 4 nGL 2 (]0,T[;F- 1 (TxM + )), 

so that in turn, 

u G C([0,r],fl'- 1 (T x E+)) nC w ([0,T];L 2 (T x M+)). 

This gives a meaning to the initial condition. 

Theorem 11.11 shows that the linearized Prandtl system (|1.5p is ill-posed in any reasonable 
sense, which strengthens the result of [4]. A difficult open problem is the extension of such 
theorem to the nonlinear setting (jl.4p . For the time being, we are unable to disprove the 
existence of a flow for the nonlinear equation. In short, we are only able to show that if the 
flow exists, it is not Lipschitz continuous from H m to H , for any m. More precisely, we 
introduce the 

Definition 1.2 (Lipschitz well-posedness for Prandtl equation near a shear flow). For any 

smooth u s and m>0, we say the IVP (jl.4p - (|1.2p is locally (H m ,H l ) Lipschitz well-posed 
if there are constants C,5o,T, and a subspace X of L°°(]0, T[; H (T x s.t. 

for any initial data uq with e y uo G H m (T x IR + ) and \\e y uo\\H m < So, there exists a unique 
distributional solution u of (|1.4j )- (|1.2j) in X , and there holds 

(1.6) ess sup |Kt)|| H i < C\\e y UQ\\ H ™ 

0<t<T ,v 

Let us quote that if u belongs to X, all terms in equation (|1.5|) are well-defined as 
distributions, including the nonlinear term 

ud x u + vd y u = d x {u 2 ) + d y (v u) G L°°(0, T; W^' P {T x R + )), Vp < 2. 

Again, the solution u satisfies 

8 t u G L°°(0,T; W~ l c ' p {T x M+)), « G C([0,T); W~J' P (T x R+)), Vp < 2, 

which gives a meaning to the initial condition. We can now state 

Theorem 1.3 (No Lipschitz continuity of the flow). There exists a shear flow u s with 
u s — U G C£°(1R+) such that: for all m > 0, the Cauchy problem (jl.4p - (jl.2p is not locally 
(ff" 1 ,^ 1 ) Lipschitz well-posed. 



1 



D. GERARD-VARET, T. NGUYEN 



The proof of the linear ill-posedness result, Theorem II, 1\ is based on the previous con- 
struction ([4]) of a strong unstable quasimode for the linearized Prandtl operator, together 
with a use of the standard closed graph theorem. For Theorem 11.31 we make a simple use 
of an idea of Guo and Tice ([6]) on deriving the ill-posedness of the flow from a strong 
linearized instability result. 

We shall end this introduction with space instability results, related to the boundary 
value problem. Let us stress that space instability is very natural in boundary layer theory. 
Indeed, many works on boundary layers are related to steady problems for flows around 
obstacles. In this context, x can be seen as the evolution variable, x = corresponding 
to the leading edge of the obstacle. Boundary layer separation, that takes place upstream 
from the leading edge, can also be seen as a blow up phenomenon in space. One can also 
refer to the paper by Metivier |10j . for a mathematical study of spacial instabilities in the 
context of Zakharov equations. 

Precisely, at the linear level, we prove the following: 

Theorem 1.4 (Spacial ill-posedness for linearized Prandtl equation). Let U > 0. There 
exists a shear flow u s with u s (y) > for y > 0, u s — U E C£°(M+) and such that : for all 
X > 0, there exists an initial data u\ satisfying 

i) e y ui G H°°(T x R + ). 

ii) The BVP (jl.4f) - (|1.3|) has no weak solution u with 

u s ueL 2 t (T;C x ([0,X];H 2 y (m + ))), u e L 2 t>x (T x (0, X); H 2 (R+)). 

This is of course an analogue of Theorem 11.11 From there, one could also obtain an 
analogue of Theorem 11.31 : we skip it for the sake of brevity. 

The proof of these spacial results are again based on construction of an unstable quasi- 
mode for the linearized operator. It turns out that the instability mechanism introduced 
in [3] to construct these unstable modes can be modified in such a way that it yields the 
ill-posedness. One should note that in course of deriving the non-existence result, we are 
obliged to prove a uniqueness result, and as it turns out, obtaining such a result is not as 
straightforward as in the case for the IVP problem. The difficulty lies in the fact that we 
now view the equation as an evolution in x and it is not at all obvious for one to obtain 
certain energy or a priori estimates for solutions of the BVP problem. Nevertheless, we 
present a proof of the uniqueness result in the last section of the paper. 

The outline of the paper is as follows: section [2] details the non-existence result stated in 
Theorem 11.11 The nonlinear result is explained in section [3j We show in section 0] how to 
adapt the arguments of [3] to get spacial ill-posedness. We shall conclude the paper with 
some comments on this spacial instability, and how it relates to the well-known results of 
Oleinik. 

2. The linearized IVP 

This section is devoted to the proof of Theorem 11.11 We start by recalling some key 
elements of article [4] . It deals with the ill-posedness of the linearized Prandtl system (jl.5p 
in the Sobolev setting. The high frequencies in x are investigated. The main point in 
the article is the construction of a strong unstable quasimode for the linearized Prandtl 
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operator C s u := u s d x u + vd y u s — d 2 u. It is achieved under the main assumption that u s 
has a non-degenerate critical point a > 0: u s (a) = 0, u"(a) < 0. More precisely, if u s takes 
the form 

\ 

(2.1) u s = u s (a) — —(y — a) 2 in the vicinity of a, 

one can build accurate approximate solutions (u™,v™) of (|1.4p that have x-frequency e^ 1 
and grow exponentially at rate e -1 / 2 . Namely, 

fj£(t,x,y) = ie i£ ~ lx e^^Ul'iy), v%(t,x,y) = e" 1 e ie ~ lx e^M*)* v?(y) 

where 

ui{e) = —u s (a) + £ 1 ' 2 t, for some r with 9mr < 0, 

and [7™, V™ are smooth functions of y. These functions are expansions of boundary layer 
type, made of 0(n) terms, with both a regular and a singular part in e. In particular, one 
has 

WKh^) > c n , \\eyU?\\ Hk{M+) < C n>fc (l + e-^ 1 )/ 4 ), VfceN 

the loss in e being due to the singular part. By accurate approximate solutions, we mean 
that d t u n e +C s u n £ = r™, with 

r?(t,x,y) = e*- 1 *^- 1 "^^, ||e^™||^ (R+) < C njfe e », Vfc. 

Again, we refer to [4] for all necessary details. Actually, the construction of [4], that deals 
with a time dependent u s , can be much simplified in the case of our steady flow u s . A 
similar construction will be described at the end of the paper, when dealing with spacial 
instability. 

We can now turn to the proof of the theorem. It ponders on the previous construction 
and on the use of the closed graph theorem. We refer to Lax [8] for similar arguments in 
the context of geometric optics. We argue by contradiction. Let us assume that for some 
T > 0, and for any data uo with e y uo E H°°(T x there is a unique solution of (jl.5p 

u e L°°(]0,r[;L 2 (T x R + )), d y u G L 2 (]0,T[xT x R + ). 

We can then define 

T : e- y H°°(T x R + ) ^ L°°(]0, T[; L 2 (T x M + )) x L 2 (]0,T[xT x R + ), u (u,d y u). 

It is a linear map between Prechet spaces, and it is easy to check that it has closed graph. 
By the closed graph theorem, we deduce that T is bounded. This means that for some 
K £N, for all u £ e -f#°°(T x M + ), 



SUp ||ti(t)||r2 + ||cLu||r2 < C 1 1 e y 1 1 . 

We quote that this bounds hold for the supremum in time and not only for the essential 
supremum, thanks to the weak continuity in time of u with values in I? . 
Let us define the family of linear operators 

S{t) : e- y H°°{T x E+) ^ L 2 (T x M+)), u (Tu ) 1 (s), te[0,T\, 
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that is S(t)uo = u(t), where u is the unique solution of (|1.5p with initial data no- The 
previous inequality shows that S(t) extends into a bounded linear operator from e~ y H K to 
L 2 , with a bound C independent of t. 

We now introduce u(t) := S(t)u™(0), and v = u — u™, where is the growing function 
defined above. It satisfies 

(2.2) dtv + C s v = -r£, v\ t=0 = 0. 

We claim that v has the Duhamel representation 

v(t) = - /** S(t-s)r*(s)ds 



Indeed, as r™ is continuous with values in e~ v H°°(T x M+), the integral is well-defined, and 
straightforward differentiation with respect to t shows that it defines another solution v of 
(12. 2D . Thus, the difference w = v — v satisfies dtw + C s w = 0, w\t=o = 0, with regularity 

w G L°°(]0,T[;L 2 (T x M + )), d y w G L 2 (0,T x T x M+). 

By our uniqueness assumption for this equation, we obtain w = 0, that is the Duhamel 
formula. 

On one hand, thanks to the bound on the S(t) and the remainder r™, we get that 
\Ht)h* < C\\eVu?(p)\\HK < Ckb- k , 

as well as 

l^, B < c£|| e ^( S )|| H K ( a )ds < C^E^He 1 ^, Vn,K>0. 
On the other hand, one has: 

|3r|t 

||«?(t)|U 2iV > c n e * 

by the properties of recalled below. Combining the last three inequalities, we deduce 

C K e- K > \Ht)\\ L 2 y > |K(t)||^ iB -|Kt)||^ iV > (c„ - C n , K s n+1 2- K ) e^r. 
This yields a contradiction for small enough e, if we take n and i such that 

,1 „ , . (ln(C ra , y /c w )+^|ln g |)^ 
2 



So far, we have proved that under assumption (|2.ip . there is for any T > an initial data 
uq G e y H°°(T x M + ) for which either existence or uniqueness of a solution of (|l,5p on [0,T] 
fails. To rule out a possible lack of uniqueness, we further assume that 

(2.3) sup ( sup | u s | + / y\d y u s \ 2 dy) < +oo. 

*>0 ^ y>0 Jo ' 

Theorem II .11 is then a consequence of the following uniqueness result: 
Proposition 2.1. Let u s be a smooth shear flow satisfying (|2.3p . Let 

w G L°°Q0, T[; L 2 (T x M+)), ^u; G L 2 (0, T x T x M + ). 
a solution of (11.51) wi/i iu|t=o = 0. Then, w = 0. 
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Proof. Let us define w k (t,y), k G Z, the Fourier transform of w(t,x,y) in x variable. We 
observe that for each k, w k solves 

d t w k + iku s w k - ikdyu s w k (y')dy' - d 2 w k = 
(2.4) I w k (t,0) = 

w k (0,y) = 0. 

Taking the standard inner product of the equation (|2.4p against the complex conjugate 
of w k and using the standard Cauchy-Schwartz inequality to the term w k dy', we obtain 

2^H^ fc lli 2 (K+) + ll 9 f* fc lli 2 (K+) - 1^1 J \ u s\\u>k\ 2 dy + \k\ J \d y u s \y 1/2 \w k \\\w k \\ L 2 {R+) dy 



< \k\( sup|it s | + / y\dyU s \ 2 dy)\\w k \ 
V t,y Jo ' 



Applying the Gronwall lemma into the last inequality yields 

\\Mt)h*(M + ) <Ce c W t \\w k (0)\\L*<jt + ), 

for some constant C. Thus, w k (t) = for each k £ Z since w k (0) = 0. That is, w = 0, and 
the theorem is proved. □ 



3. The nonlinear IVP 

With Theorem 11.11 at hand, we can turn to the nonlinear statement of Theorem 11.31 We 
shall make use of an idea of Guo and Tice [6] on deriving a nonlinear instability from a 
strong linearized one. Note that the nonlinear equation in (jl.4p reads dtu + C s u = N(u), 
with C s as in the previous section, and nonlinearity N{u) := —ud x u — vd y u. 

We now prove the theorem by contradiction. Let u s be as in Theorem 11.11 Assume that 
the IVP ([ijfD-QJl is (H m , H 1 ) locally Lipschitz well-posed for some m > 0. Let C, 8 , T be 
the constants given in the definition of Lipschitz well-posedness. By Theorem II. 1| there is 
an initial data uq G e~ y H°°(T x M + ) that does not generate any solution u of the linearized 
equation (jl.5p with 

u E L°°(]0,r[;L2(T x M+)), d y u G L 2 (]0,T[xT x M+). 

Up to multiplication by a constant, we can assume ||e y uo||_f/"' = 1. Let us take Vq := 5uq, 
with 5 a small parameter less than 5q. By the Lipschitz well-posedness hypothesis, there is 
a solution v s of (jl.4p on [0,T] with initial data Vq. Moreover, 

ess sup ||u*(t)||fl-i < CS. 
te[o,T] ,y 

In other words, u s = v s /5 is bounded in L°°(0, T; i? 1 ) uniformly with respect to 5, and 
moreover 

(3.1) d t u s + C s u s = SN(u 5 ), u s {0,x,y) =u . 

From the bound on u s , we deduce that, up to a subsequence, 

u s ^u L°° (0,^^(1 x E + )) weak * as 5 -»■ 0. 
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Furthermore, the nonlinearity 5N(u ) goes to zero strong ly in L°° (0, T; W~J' P ), for all p < 2. 
We end up with 

d t u + C s u = 0, u\ t =o = u . 

As 

u £ L°°Q0,T[;L 2 (T x R + )), 8 y u £ L 2 (]0,T[xT x R + ). 
this contradicts the result of non-existence of solutions starting from uq. 

4. Spacial instability 

This section is devoted to the boundary value problem for the linearized Prandtl equation. 
We assume U > 0, and consider some shear flow u s with u s (y) > for y > 0. As before, 
we assume u s — U £ C£°(R+), and 

ly - a) 2 

u s = u s (a) + Ug(a) in the vicinity of some a > 0, with u"(a) < 0. 

As the time t and the space x are somewhat symmetric in the Prandtl equation, we may 
adapt the construction of the unstable quasimode performed in [4]. We sketch this con- 
struction in the next paragraph, and then turn to the proof of Theorem 11.41 

4.1. The unstable quasimode. The aim of this paragraph is to construct an approximate 
solution of (|1.5p . that has high time frequency e~ l and grows exponentially for positive x 
at rate y/e . We look for growing solutions in the form 

(4.1) u(t, x, y) = e- u / £ - iul ^ x / £ u e (y), v(t, x, y) = s^e-^-^^v^y), e > 0. 

We plug the Ansatz into (|1.5p , and eliminate u £ by the divergence free condition u £ = — jj^fy . 
We end up with 



(4.2) 

Introducing the notations 



(1 + co(e)u s )v' £ — oj(e)u' s v £ — iev^ =0, y > 



Vs\ y =o = v' £]y=0 = 0. 



Cj{e) := oj{e) 1 , i := — u(e) l e, 

it reads 



(4.3) 



(u(e) + u s )v' £ - u' s v £ + ievf^ =0, y > 

V £ \y=0 = v'e\y= = 0. 



Thus, the equation gets formally close to the equation 

U 4) i + u s) v 'e - u 's v e + ievf' =0, y > 

\ Ve\y=0 = v' £ \y= = °> 

which has been studied in [3], in connection to the IVP for (|1.5p . More precisely, the authors 
build an approximate solution of (|4.4p under the form 



(4.5) 



app 



» = -u s {a) + 



v a / p (y) = H(y - a)(u s + u(e)) + e^V ( 



y-g 

eV4 
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The streamfunction Ve PP divides into two parts: a regular part 

v r / 9 (y) ■= H(y-a){u s +u{e)l 



l/2y ( vjz 



. Note that (w app (e),vr 9 ) solves (H3 



and a "shear layer part" vf{y):=e . ^ 

except for the 0{e) term coming from diffusion. However, both v r e e9 and its second derivative 
have discontinuities at y = a. This explains the introduction of vf. The shear layer profile 
V = V{z) satisfies the system 



(4.6) 



(t + <(a)y^ V' - <(a) z V + i V {3) = 0, z ^ 0, 
limV = 0. 



The jump conditions on V compensate those of the regular part. The parameter r allows 
the system not to be overdetermined. Indeed, one shows that for some appropriate r 
with 9r < 0, (14, 6h has a solution, with rapid decay to zero as z — > ±oo. This singular 
perturbation is responsible for the instability, through the eigenvalue perturbation r. We 
refer to [4j for all necessary details. Note that writing 

V(z) = V(z) + 1 K+ (V + <(a)y 

one gets rid of the jump conditions: 



(4.7) 



r + <(«)-) V 



lim V 



0. 



V 



4(a) zV + iV {3) =0, z G 
^ T + u" s (a)— as z — )• +oo 



Accordingly, the approximate solution (|4,5p takes the slightly simpler form 



(4.8) 



uj app (e) 



-u s (a) + e 1/2 r, 



i^fc/) = - a) U - «.(o) - <(a)^-^ ) + 1 \ 



y -a 



1/4 



On the basis of this former analysis, it is tempting to consider for our system (|4.3p the 
following expansion: 

-u s (a) + e 1/2 r, 

(y-ar 



(4.9) 



uffc/) = H{y -a)[u s - u s {a) - <(a) 



+ £ 



V2y 



gl/4 



However, it is not so straightforward: 

• Back to the initial notations, the first relation reads 

(4.10) uj app {e)~ l = -u s (a) + (-u app {e)~ l e) 1/2 t. 
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In particular, it is no longer an explicit definition of oo app {e). It is an equation for 
uj app (e), and we must check that the equation defines it implicitly. 
• Assuming that u app (e) is defined, e is some unreal complex number. Hence, v £ pp 
is not a priori properly defined, because V no longer has a real argument. 

Fortunately, as we will now show, these difficulties can be solved, leading to the desired 
unstable quasimode. 

We start with equation (I4,10p , A standard application of the implicit function theorem 
implies that for e small enough, the complex equation 

F(z) = 0, F(z) := z + u s (a) - (-ez) 1/2 r 

(where z 1 ' 2 is the principal value of the square root) has a unique solution near —u s (a). 
This allows to define w app (e) _1 , so co app (e). Furthermore, easy computations yield 

(4-11) co app (e) = - e 1 ' 2 —^— + o^ 2 ). 

Ug{a) u 8 (a,y/* 

Note that injecting this expression in (|4.1[) gives some exponential growth in x at rate e^ 1 ^ 2 . 

It remains to clarify the definition of v £ in (|4.9H . Of course, the complex powers e 1 / 2 , 
e 1 / 4 can still be defined. The problem lies in the profile V = V(z), which is a function over 
M, instead of (— w app (e)) _1 / 4 R. To overcome this problem, we will show that V extends 
holomorphically to the neighborhood of the real line 

U T := CX^C-r^t-oo.-llUiC-r^/^l.+oo]). 

The extension will satisfy the ODE in (j4.7j) over U T , and will have the following asymptotic 
behaviour, for small 9 > 0: 



z 2 

lim V = 0, V ~ t + <(a)— as z € e tfl R, z -»■ W e oo. 



In particular, V will satisfy all the requirements along the line (— cj(e)) 1//4 R, for e small 
enough. 

To perform our extension, we write V as 

v2- 



V{z) := [r + u:{a) Z -^W{z) 



where W satisfies 



(T + u" s {a)z 2 /2fj-W + i^(( r + <(a)z 2 /2) W) = 0, 
WmW = 0, \\m.W = 1. 

— oo +oo 



Then, performing the changes of variables 



r = ^=l<(a)| 1/2 r ; z = 2 1 /*|„J[(a)|- 1 /*5 
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leaves us with 

KmW = 0, ]haW = 1. 

— oo +00 

This system already appears in [5J. It is shown that for some appropriate f with negative 
imaginary part, there exists some solution W. Actually, much more is known on W. Indeed, 
X = W' satisfies the ODE with holomorphic coefficients 

i(f-z 2 )X" - SizX" + (t-z 2 ) 2 - 6iX = 0, 

In particular, the coefficient of the leading order term does not vanish in U T . Thus, X can 
be extended into a holomorphic solution in U T . Moreover, the previous equation on X can 
be put into the form of a first order system: 

(4.12) 4-X = zB(z)X, X = ( s -i X d x ), Biz) = ( gWW b 

In particular, B is holomorphic at infinity, with £>(oo) = (J^q). It has the two distinct 
eigenvalues ±ie l7f//4 . From there, it is possible to obtain explicit asymptotic expansions 
for X at infinity. Following [U Theorem 5.1 pl63], in any closed sector Sefi' inside which 
3f?e(ie 47r / 4 ,z 2 ) does not cancel, there are solutions X± (depending a priori on 9, 0') with the 
following asymptotic behaviour: 

(4.13) X± ~ f^^i^j e p±(i) , |z|^+oo, zeS e ,e>, 

where a± is a complex constant, and P±{z) is a polynomial of degree 2. Moreover, the 
leading term of P± is j = te " r ^ 4 ^ 2 . i n the present case, tedious computations yield 

X±{z) ~ C± z a± exp (^±ie in/A ^j , for some C + ,a+ G C, 
as \z\ — > +00, in any closed sector that is strictly inside S_n,kK 3^ , tm for some k G Z. As 

8 ~ l ~ 2 ' 8 + 2 

X decays to zero in ±00, we deduce that: for all 5 > 0, 

I J\T(5) I < C ex.p(—a\z\ 2 ), for some a > 0, z£S * , s 3n s U 5*7^ , s nn $. 

8 ' 8 8 ' 8 



This implies easily that the antiderivative W (now defined in U T ) satisfies 
and similarly 



|W^(z) — 1| < Cexp(— a'|z| 2 ), uniformly in S_ n +s 3n 

8 1 8 



< Cexp(— a'|5| 2 ), uniformly in Svn , s ug s . 

8 ' 8 

From there, one can go back to W, and then to V. 

In this way, we obtain an appropriate approximate solution of (|4.3p . Exactly as in [JJ, 
we can refine this approximation by the addition of higher order "regular" terms in Ve PP . 
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Namely, we consider an expansion of the form 

«?*(y) = H(y - a) (u s - u s (a) - <(a)^^) + e^V ^fj + £^ 2 Wv) 

Each additional term Vi >reg satisfies a first order equation of the type 

(w(e) + u s )v' i)reg - u' s v i>reg = f\ 

where f % comes from lower order terms. For example, f 2 := H(y — a)dyU s . As explained 
in [1], due to the quadratic structure of u s near a, one can show by an easy induction that 
P is identically zero near a, so that 

v i!reg = H(y - a)(u s (y) +Cu(s)) [ - J dz, 

^ ' J a {U S (Z) +0J(£)) 

is smooth across y = a. 

Back to (|4.ip . we obtain some approximate solution of (jl.5p . with frequency e^ 1 in t, 
that grows exponentially with x > at rate 

|9 : m(r)| 

a ^ ^372- 

This approximation solves p.5h up to some ^n-k^x/e 1 / 2 ^ error term in e~ y H k for all 
A;. Again, we refer to [1] for all details. 

4.2. Proof of spacial ill-posedness. Thanks to the spatially unstable quasimode from 
the last paragraph, one can try to mimic the proof of Theorem 11.11 to derive Theorem 
11.41 Indeed, all arguments related to the closed graph theorem adapt straightforwardly, as 
well as those disproving the continuity of the flow if it exists. In other words, we leave to 
the reader to check the following fact: for any X > 0, there exists some boundary data 
ui G e~ y H°°(T x M+), for which either existence or uniqueness of a solution u with 

u s ue L 2 t (T;C x ([0,X];H 2 (R + ))), u E L^{T x (0, X); 

fails. 

To complete the proof of the Theorem, it remains to show some uniqueness result for the 
BVP in this functional setting. It turns out to be more difficult than for the IVP. Therefore, 
we strengthen the assumptions on u s : besides all assumptions needed for the construction 
of the quasimode, we make the following hypothesis: the function u s , which belongs to 
U + C C °°(IR + ), satisfies 

u s (y) = y, in the vicinity of y = 0. 
Under this additional hypothesis, we are able to state 

Proposition 4.1 (Uniqueness for the BVP). Let u be a weak solution of (|1.5jl such that 

u s u e L 2 t (T;C x ([0,X];H 2 (R + ))), u € L 2 JT x (0, X); H 2 (R + )), 
and u\ x= o = 0. Then u = 0. 
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Proof. Up to a replacement of u by P^u, where P/v is the projector on temporal Fourier 
modes less than N, we can assume that 

(4.14) \d s t d^ y u{t,x,y)\ < C s \d^ y u{t,x,y)\, Vt,x,y, Va,/3. 

Our key idea to get uniqueness for equation (jl.5p is to write it as 

rv 

d x Lu = —dtu + d 2 u, where Lu := u s u — u' s I u, 

Jo 

and to use Lu and some variants of it as multipliers. Denoting by ( | ) the scalar product in 
L 2 (R + ), we have for k = 0,1, 2: 

(4-15) \d x \\d k y Lu\\ 2 L2 ^ +) = {-d t d k u\d k Lu) + (d k+2 u\d k Lu) := I k>1 + I fc>2 . 

We estimate those two terms separately, starting with = Ik,i(t, x). Let M be a constant 
such that u s = U for all y>M. We claim that : Ve > 0, V/c = 0, 1, 2, 

(4.16) J h,x{t)dt<C £ (\Ht)\\l HQM) + \\d y u{t)\\l^ M) ) dt+e \\d*u(t)\\l 2{0>M) dt. 

Note that, for clarity, we omit to indicate the dependance on x of the various quantities. 
For brevity, we just prove our claim in the case k = 2, which is the most involved: one has 

/2,i = — (dtd y u\d y Lu) = — I dtdyU I u s dyU + u s dyU — uy' / u — u" s u ) Ay 

Jr+ \ Jo J 

= ~\~u I u s \d y u s \ 2 dy- I dt&u(u' a d y u-uf) [ u-u" s u]dy. 
1 Ut J R+ Jr + v Jo J 

Integration in time makes the first term at the r.h.s. vanishes. The other terms involve 
derivatives of u s , so that the integrals over M + can be replaced by integrals over (0,M). 
Using Cauchy-Schwarz inequality, we end up with 

h,\ < C / ||9t^n|| L 2 (0iM ) I 11^14^2(0^/) + || / -u|| L 2( 0)M) + ||u|| £ ,2( 0jAf) 



< C / \\dtd 2 u\\ L 2 {QM) (\\d y u\\ L 2^ M) + (l + M 2 /2)|M| L 2 {0jM) ) 

< e / \\ d y u \\h(o,M) + ° £ / (W d y u \\LHo,M) + \H\l^{o,m)) 2 

where Young's inequality and (|4.14p were used for the last inequality. This proves the claim. 

We now turn to the estimate of I k %■ Let m > such that u s (y) = y for y < m. We claim 
that: for all k = 0, 1,2, 

(4.17) I fcj2 < C (||w||| 2(0iM) + \\dyu\\ 2 L2(0M) + ||aJu||ia (m)Af) ) • 

Again, we omitted the t, x dependence in the notations. These inequalities follow from 
simple integration by parts. As before, we consider the case k = 2, the other ones being 
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easiest. We have 



y u — 



u — u" s u I dy 



h,2 = {d 4 y u\d 2 y Lu) = I tfu (u s d 2 u + u'A 

Jr+ V 

= - / dlud y (u s d 2 y u + u' s d y u)dy - [ d 2 y ud 2 y (u^ P ' u + u" s u]dy 
Jr+ Jr+ V Jo J 

through integrations by parts. The last term at the r.h.s involves derivatives of u' s , so is for 
m < y < M . It also involves up to second order derivatives of u only. Thus 

h,2 < - d^udy (u s d 2 u + u' s d y u) dy + C^\\d^u\\ 2 L2{mM) . 
+ k=0 
As regards the first term, we have 

dyudy (u s d 2 u + u' s d y u) dy = — dyU (u s dyU + 2u' s d 2 u + u"d y u) 

Jr + 

d z y u (2u' s d 2 u + u" s d y u) = - j u' s dy(d 2 u) 2 + f d 2 y ud y {u% 
r 2 

4(d 2 u) 2 + / d 2 y ud y {u" s d y u) < C^||^u||| 2 
Jr+ k=1 



-I 



< - 



U 



J 

JR- 



2 (m,M)- 



This proves our claim. 

We must now link the Sobolev norms of u and Lu. This is the purpose of the 
Lemma 4.2. Let u = u(y) G H 2 (0,M), u\ y=0 = 0. Then, 

IMI/^tW) < C \\ Lu \\h 2 (o,m) 



and 



\ u \\H 2 (m,M) 



< C m \\Lu\ 



H 2 (0,M) 



Proof of the lemma. Note that by our assumptions on u s , Lu belongs to H 2 (0,M), and 
satisfies Lu\ y= o = d y Lu\ y= o = 0. Denoting / := Lu, we have an explicit representation of u 
in terms of /, by solving the first order ODE Lu = f: 



<v) = «.(y) / —dt + ^ 



4(v) / 

Jo 



« S (*) 

78 /(*) 



t 2 



' /•" /(*) ,,/ . /(f) 



+ u 2 (y) + u 3 (y). 



Clearly, \\u 2 + u 3 \\ H 2( mM) < C m ||/|| H 2( 0) M)- As regards ui, we compute 



Jo * Jo 



to 



so that 

II u i||h 2 (o,m) < c 



/' 

JO 



i 2 



dt 



< c 



fit) 



1/2 



+ 



L°°(0,M) 



< C" 



H 2 (0,M)> 
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using Hardy and Sobolev inequalities. We thus obtain the second inequality. For the proof 
of the first inequality, it remains to control u on the interval (0, m). It satisfies there the 
equation yu — f Q y u = f , which gives yu' = /', and from Hardy's inequality, we deduce that 
ll u/ |lL 2 (o,m) — ll/ll_ff 2 (o,m)- This ends the proof of the lemma. 

Combining this lemma with (|4.15p . (|4.16p and (|4.17|) . we obtain that the solution u of 
(|1.5p with boundary data u\ x =o = satisfies 



(4.18) \\Lu(x)\\ 2 L 2 {T . H 2 {R+)) < C £ J \\Lu{x')\\ 2 L 2 {T . H ^ R+)) dx' + £||d>||i 2((0>x)xTxR+) 

for all e > (e = 1 will be enough here). We still need to control the 1? norm of d 2 u. 
Therefore, we differentiate (jl.5p with respect to y, and write the resulting equation as 



u s d x L'u = d^u — dtd y u, L'u := d y u — — / u. 



Us Jo 

Note that the second term defining L'u has no singularity at zero, as u" vanishes identically 
for y < m. Multiplying this equation by L'u, and integrating in t, x, y we obtain 

/ it \ I 1 /l i r 1 I ~ — / ( 1 r ^ , 



2 

One has easily, still with (|4.14|) 



u s \L'u{x)\ 2 = / (dyU — dtd y u)L'u. 

TxR + J(0,x)xTxI 



/ {-d t d y u)L'u < cj2\\d k y n\\h mx 



'(0,i)xTxt + k=Q 

whereas a simple integration by parts yield 

2 



)xTx(0,M))' 



/ d^uL'u = - \d 2 u\ 2 + ^||^«||| 2(((! X! ,„ , /n 

i(0,rr)xTxR + J(0,x)xTxR + k=Q 



Hence, 

|2j„,j + , ||»2 



(4.19) - / u s \L'u{x)\ 2 dydt + \\d 2 u\\ L 2 {{Q 

,i)x)TxR + ) 

z ./TxR + 



- C I 7, \\®U U \\l 2 ((0.x)xTx(0.M)) + II^Hli 2 ((0,a;)xTx(m,.M")) 
\k=0 

Combining this inequality with Lemma 14.21 implies that 

ll^ U lli 2 ((0,x)xTxR + ) - < ^W Lu \\l 2 ((Q,x)xT;H2(R + )) 

Together with (|4.18j) . this leads to 

2 ^ r< I II r_./_'\il2 



\\ Lu {x)\\ H 2 y( j xR+) < C e J \\Lu(x')\\ L 2( T . H 2 {R+) )dx' . 

By Gronwall lemma, we deduce that Lu = 0, and from there that u = 0. This concludes 
the proof of the proposition, and in turn the proof of Theorem ll.4i □ 
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4.3. Final comments. The ill-posedness results discussed in this paper may give some 
insight into the classical results of Oleinik on boundary layer theory. These results, collected 
in the book [11], establish various well-posedness theorems for the Prandtl equation, under 
some monotonicity properties for the data with respect to y. Interestingly, the monotonicity 
assumptions used by Oleinik are different whether the steady case or the unsteady one is 
considered. Roughly: 

• In the steady case, the main hypothesis is that u| x =o be monotonic near 0, that is 
d y u\ x=0 (0) > 0. 

• But in the unsteady case, the monotonicity hypothesis is stronger. The data u\ x= o 
and u\t=o are assumed to be monotonic everywhere, that is d y u\t=o(x, y) > and 
dyu\ x =o(t,y) > for all y > 0. 

In Oleinik's work, these assumptions appear as technical, connected to some special changes 
of variables that are different in the steady and unsteady case (these are the so-called Von 
Mises and Crocco transforms). 

From our results, one may believe that this difference in the steady and unsteady as- 
sumptions is more than technical. Indeed, we have proved that as soon as time is involved 
in the Prandtl system, either through an IVP or a BVP with time periodicity, a lack of 
monotonicity somewhere in the flow (u' s (a) > for some a > 0) is enough to trigger a strong 
instability (at least at the linearized level). On the contrary, as soon as a steady situation 
is considered, this instability mechanism dissapears. For instance, the space instability de- 
velopped in this section relies on high frequency time oscillations, and is not there when 
the steady BVP is considered. Actually, the BVP for the steady version of (|1 . 5[) : 



can be shown to be well-posed for the shear flow u s considered in Proposition 14.11 Indeed, 
the estimates established there to prove uniqueness can be used in the same way to establish 
some a priori estimates, and from there a well-posedness result (in appropriate functional 
spaces). Again, this is coherent with the works of Oleinik, as we assumed in Proposition 14. II 
that u s (y) = y for small y. Still in the spirit of this proposition, it would be very interesting 
to recover the well-posedness results of Oleinik without special changes of variables, that is 
simply through linearization and/or energy estimates. 
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